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Abstract 

The Schrodinger problem of deducing the microscopic dynamics from the input- 
output statistics data is known to admit a solution in terms of Markov diffusions. The 
uniqueness of solution is found linked to the natural boundaries respected by the 
underlying random motion. By choosing a reference Smoluchowski diffusion, we au- 
tomatically fix the Feynman-Kac potential and the field of local accelerations it 
induces. We generate the family of affiliated diffusions with the same local dy- 
namics, but different inaccessible boundaries on finite, semi-infinite and infinite do- 
mains. For each diffusion process a unique Feynman-Kac kernel is obtained by the 
constrained (Dirichlet boundary data) Wiener path integration. As a by-product of 
the discussion, we give an overview of the problem of inaccessible boundaries for 
the diffusion, and bring together (sometimes viewed from unexpected angles) results 
which are little known, and dispersed in publications from scarcely communicating 
areas of mathematics and physics 
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1 The Schrodinger problem: microscopic dynamics from 
the input-output statistics 

According to M.Kac' 1 ' any kind of time developement (be it deterministic or essentially 
probabilistic), which is analyzable in terms of probabilitydeserves the name of the stochas- 
tic process. 

Given a dynamical law of motion (for a particle as example), in many cases one 
can associate with it (compute or approximate the observed frequency data) a probabil- 
ity distribution and various mean values. In fact, it is well known that inequivalent finite 
difference random motion problems may give rise to the same continuous approximant 
(e.g. the diffusion equation representation of discrete processes). As well, in the study of 
nonlinear dynamical systems' 2 ' ,given almost any (for the purposes of our discussion , ba- 
sically one-dimensional)probability densityit is possible to construct an infinite number 
of deterministic finite difference equations,whose iterates are chaotic and which give rise 
to this a priori prescribed density. 

The inverse operation of deducing the detailed (possibly individual, microscopic) 
dynamics, which either implies or is consistent with the given probability distribution 
(and eventually with its own time evolution) cannot thus have a unique solution. 

If we disregard the detailed nature (like its chaotic Jump pro cess, random walk,phase 
space process with friction etc. implementations) of the given process, it appears' 3,4 ' that 
the standard Brownian motion and/or the broad class of Markovian diffusions incorpo- 
rating the Wiener noise input, provide satisfactory approximations for a large variety of 
phenomena.lt especially pertains to the explicit modelling of any unknown in detail phys- 
ical process in terms of the input-output statistics (conditional probabilities and averages) 
of random motions with a finite time of duration. 

From now on ,we shall confine our attention to continuous Markov processes, whose 
random variable X(t),t > takes values on the real line i? 1 ,and in particular can be 
restricted (constrained) to remain within the interval A C R 1 , which may be finite or 
(semi-) infinite but basically an open set. Its boundaries <9A (endpoints) will be denoted 
ri,r 2 with — oo < r\ < r 2 < oo. 

In the above input-output statistics context, let us invoke a probabilistic problem,ori 
ginally due to Schrodinger' 5-7 ] : given two strictly positive (on an open interval) boundary 
probability distributions p (x), px(x) for a process with the time of duration T > O.Can 
we uniquely identify the stochastic process interpolating between them ? 

Perhaps unexpectedly in the light of our previous comments ,the answer is known' 6 ' 7 ' 
to be affirmative ,if we assume the interpolating process to be Markovian. In fact,we get 
here a unique Markovian diffusion,which is specified by the joint probability distribution 
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J dym(x,y) = p (x) 
dxm(x,y) = p T {y) 



where 

m(x, y) = 0* (x, 0) k(x, 0, y, T) G(y, T) (2) 

and the two unknown (not necessarily Lebesgue integrable) functions Q*(x,0),Q(y,T) 
come out as solutions of the same sign of the integral identities (1). Provided ,we have at 
our disposal a bounded strictly positive integral kernel k(x, s, y, t), < s < t < T.Then: 

e,(x, t) = J fc(0, y, x, *)6,(y, 0)rfy (3) 

6(x, s) = J k(s,x,y,T)e(y,T)dy 

and the sought for interpolation has a probability distribution p(x,t) = (6*0) (x, t), t G 
[0,T].The transition density 

p(y, s, x, t) = k(y, s, x, t) ® X ' t (4) 

Q(y, s) 

with s < t is a fundamental solution (p — > 5(x — y) as t | s) of the forward Kolmogorov 
(i.e. Fokker-Planck) equation with a diffusion constant D > (this choice narrows slightly 
the allowed framework): 

d t p = DA x p-V x (bp) (5) 



p(x, t) = J p(y, s, x, t) p(y, s) dy 
with po(x) = p(x,0) and the drift b(x,t) given by : 

V0 

b(x,t) = 2D—(x,t) (6) 
The backward diffusion equation is solved by the same transition density 

d s p = —DAyp - bV y p (7) 

p = p(y,s,x,t) , s<t, b = b(y,s) 

and we deal here with a unique diffusion process, whose transition density is a common 
fundamental solution for both the backward and forward Kolmogorov equations. 
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To understand the role of the integral kernel k(y, s, x, t) in (l)-(7) let us assume that 
Q(x,t) is given in the form (drifts are gradient fields as a consequence): 

Q(x, t) = ±exp$(x, t) b(x, t) = 2DV$(x, t) (9) 

x e (ri,r 2 ) 

and insert (4) to the Fokker-Planck equation (5).Then[ 8 ' 13 ' 14 l, if p(y,s,x,t) is to solve 
(5), the kernel k(y,s,x,t) must be a fundamental solution of the generalised diffusion 
equation: 

dt k = DA X k - — Sl(x, t) k (10) 
k(y, s, x, t) — > S(x — y) ast { s 

n(x,t)=2D[d t $ + l^ + Vb)] 

and to guarrantee (3) ,it must display the semigroup composition properties. 

Notice that (4), (9) imply that the backward diffusion equation (7) takes the form 
of the adjoint to (10): 

d s k = -DA y k + ^n(y,s)k (11) 
k = k(y, s, x, t) 

If the process takes place in-between boundaries at infinity r\ = — oo, r 2 = +oo ,the 
standard restrictions on the auxiliary potential CI (Rellich class' 15 ' 16 ]) and hence on the 
drift potential $(x, t) , yield the familiar Feynman-Kac representation of the fundamental 
solution k(y,s,x,t) common for (10) and (11): 

k(y, s, x, t) = J exp[-7^j J g Ct(X (u) , u)du] dfi[s, y\t,x] (12) 

which integrates exp[— (1/2D) /* Q(X (u) , u)du] weighting factors with respect to the con- 
ditional Wiener measure i.e. along all sample paths of the Wiener process which connect 
y with x in time t — s.See e.g. also Refs. 17,18. More elaborate discussion is necessary ,if 
at least one of the boundary points is not at infinity. 

Let us notice that the time independence of fl is granted if either $ is independent 
of time ,or depends on time at most linearly. Then the standard expression exp[—H(t — 
s)](y,x) for the kernel k clearly reveals the involved semigroup properties,with H = 
—DA + (1/2D)Vl(x) being the essentially self adjoint operator on its (Hilbert space) 
domain. 
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2 Natural boundaries make diffusions unique 



We shall make one more step narrowing slightly the scope of our discussion by admitting 
diffusions (l)-(7) whose drift fields are time-independent:<9t&(x, t) = for all x.We know' 8 ' 
that both the free Brownian motion and the Brownian motion in a field of force in the 
Smoluchowski approximation , belong to this class of processes. We know also' 9 ' that the 
boundary value problems for the Smoluchowski equation have a profound physical signifi- 
cance, albeit the attention paid to various cases is definitely unbalanced in the literature.lt 
is then interesting to observe that the situation we encounter in connection with (l)-(7) 
is very specific from the point of view of Feller's' 9-12 ' classification of one-dimensional 
diffusions encompassing effects of the boundary data. Our case is precisely the Feller dif- 
fusion respecting (confined between) the natural boundaries .An equivalent statement is 
that boundary points n, r 2 are inaccessible barriers for the process i.e. there is no positive 
probability that any of them can be reached from the interior of (ri,r 2 ) within a finite 
time for all X(0) = x e (ri,r2),see e.g. Refs.10,11 and 12, chap. III. 4. 

In the mathematical literature' 10,11,18 ' a clear distinction is made between the back- 
ward and forward Kolmogorov equations. The backward one defines the transition density 
of the process,while the forward (Fokker-Planck) one determines the probability distribu- 
tion (density) of diffusion. With a given backward equation one can usually associate the 
whole family of forward (Fokker-Planck) equations, whose explicit form reflects the par- 
ticular choice of boundary data.This fundamental distinction seemingly evaporates in our 
previous discussion (l)-(ll) ,but it is by no means incidental. In fact , according to Feller' 10 ': 
in order that there exists one and only one (homogeneous ; p(y,s,x,t) = p(t — s;y,x)) 
process satisfying — dtu = DAu + U\7u in a finite or infinite interval r\ < x < r 2 it is 
necessary and sufficient that both boundaries are inaccessible (the probability to reach 
either of them within a finite time interval must be zero). 

A general feature of the inaccessible boundary problems is that the density of diffu- 
sion vanishes' 9 - 20 ' at the boundaries:p(rx) = = p(r 2 ). This property is shared with ,more 
familiar in the realm of the statistical physics, absorbing barrier processes. The link is in- 
deed very close' 10 ': conventionally the absorbing barrier process is defined on the closed 
interval [Ri,R 2 ] , however we can always consider it on the open set (i?i, R 2 ) .Theorem 7 
of Ref.10 states that ,if the boundaries n < R\ < R 2 < r 2 are inaccessible for the process 
,then transition densities of the absorbing barrier process on (i?i, R 2 ) as R\ — > n, R 2 — > r 2 
converge to the unique (!) transition density of the diffusion with unattainable boundaries 
on (ri,r 2 ). It implies that locally ,the inaccessible boundary problem , in principle can 
be modelled (approximated) to an arbitrary degree of accuracy by the absorbing barrier 
process. 

It is interesting to notice that the classification of Feller's boundaries in the homoge- 
neous case (time independent drifts, continuous but not necessarily bounded in the interval 



5 



(fi,f2)) follows from investigating the Lebesgue integrability of the Hille functions (see 
Refs.9,11,20 for more details) : 



where x G (ri,r 2 ) and we have used (9) in the above. Apparently ,the <&(xq) contribution 
in Li(x) is irrelevant, and the integrability of exp[±<3>(x)] matters here. 

If Li(x) is not Lebesgue integrable on [x , R], where R = r\ or r 2 ,then R stands for 
the natural repulsive boundary of the diffusion. If Li(x) is integrable but L 2 (x) is not ,then 
R is a natural attractive boundary for the process. both being inaccesible .As indicated in 
Ref.9 there is no universally established terminology and a certain discrimination between 
Feller's and the Gihman-Skorohod definition is possible, albeit without consequences for 
our discussion. 

Following Ref.12 let us denote P x [tr < oo] a probability that a process originating 
from x G (ri,r 2 ) would hit the point R at the moment tr for the first time. Then, the 
inaccessibility of boundaries can be expressed by the statements: 



For the natural boundary ,R is inaccessible from the interior of (r 1? r 2 ) and the interior of 
(r 1 ,r 2 ) is inaccessible from i?.Following this terminology: (i) An inaccessible boundary is 
called attractive ,if for any e > there exists 5 > such that P :r [/im^ 00 X(t) = R] > 1 — e 
for all x G (R, R + 5) in case of the left , while x G (R — S, R) in case of the right boundary, 
(ii) An inaccessible (left) boundary is called repelling ,if for any x > R and y < x we have 
P y [r x < oo] = 1. 

Remark 1: The standard (unrestricted ) Brownian motion on R 1 is the most 
obvious example of diffusion with natural boundaries.lt is not quite trival to construct 
explicit examples ,if one of the boundaries is not at infinity. The classic example of diffusion 
on the half-line with natural boundaries at and +oo is provided by the so called Bessel 
process^ 11 ' 21 ] ,with the diffusion (backward) generator L a = A x + V x (we absorb the 
diffusion constant D in the rescaled time parameter). The point r\ = is never reached 
with probability one if a > 0. In case of a = ,the transition density reads' 21 ': 




(13) 




right boundary, Vx < R , P x [tr < oo] = 



left boundary , \fx > R , P x [tr < oo] = 



pit; Xq, x) = const— exp[ 



x 




(14) 
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where the modified Bessel function (Ref.22,chap.7) is given by Io(a) = Yl'jLo Q;2j /[2 2 ' 7 O0 2 ]- 
Another (less explicit, as given in the form of estimates for the transition density) example 
pertains' 23 ^ to the diffusion equation with the one dimensional harmonic oscillator poten- 
tial on the half-line x > a > .The related constrained path integrals are considered in 
Refs.24,25. 

As mentioned before , diffusions with inaccessible barriers might have drifts which 
are unbounded on (r±, r 2 ). Hence ,our discussion definitely falls into the framework of 
diffusion processes with sing ular drift fields [26 ' 27] ,see also Refs. 28-33, which is not covered 
by standard monographs on stochastic processes. Particularly illuminating in this respect 
is the analysis of Ref.26 where for quite general diffusions, the unattainability of nodal 
sets (on which the probability density vanishes) in a finite time was demonstrated ,in 
the sense that P x [tr = 00} = 1. This crucial property (valid for diffusions with natural 
boundaries as well) allows to extend the theory of stochastic differential equations and 
integrals to diffusions , whose drifts show up a bad (unboundedness or divergence to 
infinity) behaviour when approaching the boundaries. 

We skip the standard details concerning the probability space , filtration and the 
process adapted to this filtration (see however Refs. 22, 34, 35) and notice that a continuous 
random process X(t),t G [0, T\ with a probability measure P is called a process of the 
diffusion type if its drift b(x) obeys: 

P[ [ T I b(X(t)) I dt < 00] = 1 (15) 
Jo 

and ,given the standard Wiener process (Brownian motion) W(t) ,the integral identity (D 
constant and positive) holds true P-almost surely (except possibly on sets of P-measure 
zero). It means that W(t) = (l/v / 2~D[X(t) — J* * b(X(s))ds] is a standard Wiener process 
with respect to the probability measure P of the process X(t). 

For diffusions with natural boundaries , we remain within the regularity interval of 
b(X(t)) for all (finite) times ,and (15) apparently is valid. Therefore, the standard rules 
of the stochastic Ito calculus' 20 ' can be adopted to relate the Fokker-Planck equation (7) 
with the natural boundaries to the diffusion process X(t), which' 35 ] "admits the stochastic 
differential" 

dX(t) = b(X(t))dt + V2DdW(t) (17) 
X(0) = x o ,te[0,T\ 

for all (finite ) times. The weak (in view of assigning the density po(x) to the random 
variable A(0)) solution of (17) is thus well defined. 

For stochastic differential equations of the form (17) ,the explicit Wiener noise input 
,because of (9) implies that irrespective of whether natural boundaries are at infinity or 
not ,the Cameron-Martin-Girsanov' 37-39 ! method of measure substitutions which parallel 
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transformations of drifts ,is applicable. Even though the drifts are generally unbounded 
on (ri,r 2 ) and the original theory' 37-39 ^ is essentially based on the boundedness demand. 
It is basically due to the fact that the probabilistic Cameron-Martin formula relating 
the probability measure Px of X(t) with the Wiener meaasure Pyy (strictly speaking 
it is the Radon-Nikodym derivative of one measure with respect to another) reduces to 
the familiar Feynman-Kac formula' 17 ' 18 ' 8,14,30,36 ! (with the multiplicative normalisation). 
The problem of the existence of the Radon-Nikodym derivative (and this of the absolute 
continuity of Px with respect to Pw , which implies that sets of Piy-measure zero are of 
Px-measure zero as well) is then replaced by the standard functional analytic problem' 15 ' 16 ' 
of representing the semigroup operator kernel via the Feynman-Kac integral with respect 
to the conditional Wiener measure. 

The Feynman-Kac formula is casually viewed to encompass the unrestricted (the 
whole of R n ) motions,however it is known to be localizable ,and its validity extends also 
to finite and semi-infinite subsets of R 1 (R n more generally) as demonstrated in the context 
of the statistical mechanics of continuous quantum systems' 16 ' 23 ' 40_44 '.More specifically, it 
refers to the Dirichlet boundary conditions for self-adjoint Hamiltonians, which ensure 
their essential self-adjointness (to yield the Trotter formula ). 

Remark 2 : It is perhaps worthwhile to say few words about the situation when the 
process , in principle can reach or cross the boundary (nodal surface of the probability 
density in higher dimensions) in a finite time. As before we limit our discussion to the 
stationary Markov diffusion process and assume that the drift b is a gradient ~ Vp/p.If 
p = exp$ where p 1 / 2 is not an element of the Sobolev space H} oc (R d ) (i.e. is not integrable 
on bounded sets, with its first derivative), then it is known ' 49 1 that the process can reach or 
cross the set N+ = [x E R d | p = 0] It is possible' 63 ] to formulate a useful criterion (easy 
to generalize to higher dimensions) for a "tunelling" (transmission) through a chosen 
point in i? 4 ,let it be the origin x — 0. Let us take b ~ V@/O.The dynamics is given 
by the energy form on L 2 (R,dx) i.e. E[f±, f 2 ) = J"V fiV f2p(x)dx where obviously p = 
00* = 2 , compare e.g. (l)-(7). If p{x) < A \ x | in the close vicinity of the origin 
,at least on one side ,then there is no tunelling through x = .The process cannot 
cross this point ,but may be absorbed or reflected . One needs a bit stronger restriction 
to prevent the diffusion from hitting the node : non-transmission and non-hitting are not 
equivalent concepts, although there is no communication between diffusions on the positive 
and negative semi-axis respectively. On the other hand ,if p{x) > A\ x \ a with < a < 1 
in a neighbourhood of the origin, then there is a particle transmission (tunelling) through 
x = O.As observed in Ref.63 , according to Feller's classification of boundaries, stands for 
the regular boundary in this case. 
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3 Hydrodynamic representation:local conservation laws 
and the Newtonian dynamics in mean 

Let us emphasize the importance of (17), and of the Ito differential formula induced by 
(17) for smooth functions of the random variable X(t). Its first consequence is that given 
p(y,s,x,t), for any smooth function of the random variable the forward time derivative 
in the conditonal mean can be introduced^ 4 ' 7,20 ' 27 ' (we bypass in this way the inherent 
non-differentiability of sample paths of the process) 

Hm Atl0 — [ J p(x, t,y,t + At)f(y, t + At)dy - f(x, t)} = (D + f(X(t),t) = (18) 

= (dt + bV + DA)f(X(t),t) 
X(t) = x 

so that the second forward derivative associates with our diffusion the local field of accel- 
erations: 

(D 2 + X)(t) = (D + b)(X(t),t) = (dtb + bVb + DAb)(X(t),t) = Vn(X(t),t) (19) 

with the (auxiliary potential Q(x,t) introduced before in the formula (10). Since we have 
given p(x,t) for all t G [0, T], the notion of the backward transition density p*(y, s, x, t) 
can be introduced as well 

p(x, t)p*(y, s, x, t) = p(y, s, x, t)p(y, s) (20) 

which allows to define the backward derivative of the process in the conditional mean 
(cf.Refs.8,45-47 for a discussion of these concepts in case of the most traditional Brownian 
motion) 

lim Atl0 ^ t [x- J p*{y,t-At,x,t)ydy] = (D-X)(t) = K(X(t),t) = [b-2DVlnp](X(t),t) 

(21) 

(D_f)(X(t),t) = (d t + K V - DA)f(X(t),t) 

Apparently ,the validity of (19) (cf.Refs.7,14,27 for related considerations) extends 
to (D 2 _X){t) as well ,and there holds 

(D 2 + X){t) = (D 2 _X)(t) = d t v + vVv + VQ = Vn (22) 
v(x,t) = ]-(b + b*)(x,t) , u(x,t) = ^(b - b*)(x,t) = DVlnp(x,t) 
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Q(x,t) = 2I?-± r 

Clearly ,if b and p are time-independent ,then (22) reduces to the identity 

vVv = V(Q - Q) (23) 

while in case of constant (or vanishing) current velocity v, the acceleration formula (22) 
reduces to 

= V(fi - Q) (24) 

which establishes a very restrictive relationship! 48-52 ] between the auxiliary potential Q(x) 
(and hence the drift b(x)) and the probability distribution p(x) of the stationary dif- 
fusion. The pertinent random motions have their place in the mathematically oriented 
literature' 48-53 ]. 

Let us notice that (22) allows to transform the Fokker-Planck equation (7) into 
the familiar continuity equation, so that the diffusion process X(t) admits a recasting in 
terms of the manifestly hydrodynamical local conservation laws (we adopt here the kinetic 
theory lore) 

d tP = -V(pv) (25) 

d t v + vVv = V(ft - Q) 

Po(x) = p(x, 0) , v (x) = v(x, 0) 

which form a closed ( in fact,Cauchy) nonlinearly coupled system of differential equa- 
tions, strictly equivalent to the previous (7), (19). 

In view of the natural boundaries (where the density p(x, t) vanishes) ,the diffusion 
respects a specific (" Euclidean looking" ) version of the Ehrenfest theorem! 14 ] : 

E[VQ] = (26) 

^E[X(t)] = jE[v{X{t),t)\ = E[(d t v + vVv)(X(t,t)} = E[Vn(X(t),t)} 

Notice that the auxiliary potential of the form Q = 2Q — V where V is any Rellich class 
representative, defines drifts of Nelson's diffusions^ 14 ' 47 ] for which E[VQ] = =>■ E'fVfi] = 
— E'fVV] i.e. the "standard looking" form of the second Newton law in the mean arises. 

At this point it seems instructive to comment on the essentially hydrodynamical 
features (compressible fluid/gas case)of the problem (25), where the "pressure" term VQ 
is quite annoying from the traditional kinetic theory perspective' 54 ' 55 '. Although (25) has a 
conspicuous Euler form ,one should notice that if the starting point of our discussion would 
be a typical Smoluchowski diffusion' 8 ] (7), (17) whose drift is given by the Stokes formula 
(i.e. is proportional to the external force F = — W acting on diffusing molecules), then its 
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external force factor is precisely the one retained from the original Kramers phase-space 
formulation^ 3 ' 4,9 ] of the high friction affected random motion. In the Euler description of 
fluids and gases , the very same force which is present in the Kramers (or Boltzmann in 
the traditional discussion) equation , should reappear on the right-hand-side of the local 
conservation law (momentum balance formula) (25) .Except for the harmonic oscillator 
example ,in view of (10) it is generally not the case in application to diffusion processes. 

Following the hydrodynamic tradition let us analyze the issue in more detail. We 
consider a reference volume (control interval) [a,/3] in R 1 (or A C R 1 ), which at time 
t G [0, T] comprises a certain fraction of particles (fluid constituents) ,for an instant of 
course. Since we might deal with a flow (proportional to the current velocity v(x,t),(22)) 
the time rate of particles loss by the volume [a, (5] at time t ,is equal to the flow outgoing 
through the boundaries i.e. 

rP 

-d t p(x,t)dx = p(0,t)v(P,t) - p(a,t)v(a,t) (27) 

J a 

which is a consequence of the continuity equation. To analyze the momentum balance, let 
us allow^ for an infinitesimal deformation of the boundaries of [a,f3] to have entirely 
compensated the mass (particle) loss (27) 

[a, /?]->[« + v(a, t)At, (3 + v((3, t) At] 

Effectively, we pass then to the locally co-moving frame .It implies 

UrriAtio — [ p(x,t + At)dx — p(x,t)dx] — (28) 

J a+v a At J a 



A, L / 
t Ja+v a At 

— UmAtw-r- [ / p(x,t)dx + At / (d t p)dx + / p(x, t)dx] = 

At Ja+v a At Ja J/3 

Let us investigate what happens to the local flows (pv) (x, t) , if we proceed in the same 
way (leading terms only are retained): 

r l3+v At r /3 

/ (pv)(x,t + At)dx - (pv)(x,t)dt~ (29) 

Ja+v a At Ja 

~ -(pv 2 )(a,t)At + (pv 2 )((3,t)At + At / [d t (pv)}dx 
Because of (25) we have 

dt(pv) = -V{pv 2 ) + p\7(tt - Q) (30) 
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and the rate of change of momentum associated with the control volume [a,/3] is 

I r/3+VpAt r p r p 

Um Atl0 — [ (pv)(x,t + At) - / (pv){x,t)} = P V(Sl-Q)dx (31) 

lAt Ja+VaAt Ja J a 

However' 45 ^ 

VP 

VQ = (32) 

P 

P = D 2 pAlnp 

and consequently 

/ p V(Q - Q)dx = f pVVtdx - f VPdx = (33) 

J a J a J a 

= E[VQf a + P(a,t)-P(f3,t) 

Clearly , Vil refers to the Euler-type volume force , while VQ (or more correctly , P) refers 
to the "pressure" effects entirely due to the particle transfer rate through the bound- 
aries of the considered volume. The latter property can be consistently attributted to the 
Wiener noise proper:it sends particles away from the areas of larger concentration. See 
e.g. Refs. 14,45-47 and also Ref.46 for a discussion of the Brownian recoil principle, which 
reverses the originally Wiener flows. 

As it appears ,the validity of the stochastic differential representation (17) of the dif- 
fusion (5) implies the validity of the hydrodynamical representation (25) of the process.lt 
in turn gives a distinguished status to the auxiliary potential Q(x,t) of (10)-(12).We 
encounter here' 8,14 ' a fundamental problem of what is to be interpreted by a physicist 
(observer) as the external force field manifestation in the diffusion process. Should it be 
dictated by the drift forml 3,4,9,57 ' following Smoluchowski and Kramers , or rather by Vfi 
entering the evident (albeit "Euclidean looking") second Newton law , respected by the 
diffusion ? In the standard derivations of the Smoluchowski equation ,the deterministic 
part (force and friction terms) of the Langevin equation is postulated. What however, if 
the experimental data pertain to the local conservation laws like (25) and (27), and there 
is no direct (experimental) access to the microscopic dynamics ? 

If the field of accelerations Vf2 is taken as the primary defining characteristics of 
diffusion we deal with ,then we face the problem of deducing all drifts ,and hence diffu- 
sions , which give rise to the same acceleration field, and thus form a class of dynamically 
equivalent diffusions . 
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4 Feynman-Kac kernels on finite and semi-infinite 
domains with inaccessible boundaries: dynamically 
equivalent diffusions 

Let us analyze the second consequence of the unattainability of the boundaries , which via 
(15) gives rise to (17). On the same footing as in case of (15) ,we have satisfied another 
probabilistic identity: 

P[ [ b 2 (X(t))dt < oo] = 1 (34) 
Jo 

For a diffusion X(t) with the differential (17) , Theorem 6 of Ref.35 states that (34) is 
a sufficient and necessary condition for the absolute continuity of the measure P = Px 
with respect to the Wiener measure P w . Since ,for any (Borel) set A ,P W (A) = implies 
Px{A) = , the Radon-Nikodym theorem applies' 34 ' and densities of these measures can 
be related.lt is worthwhile to mention the demonstration due to Fukushima' 49 ' that the 
mutual absolute continuity (the previous implication can be reversed) holds true for most 
measures we are interested in. 

In the notation (12), the conditional Wiener measure d/i[s,y \ t,x] gives rise to the 
familiar heat kernel, if we set Q = identically.lt in turn induces the Wiener measure Pw 
of the set of all sample paths, which originate from y at time s and terminate (can be 
located) in the Borel set A after time t — s: 

Pw[A] = J dx J d(j,[s,y \t,x\ — J d/i (35) 

where ,for simplicity of notations, the (y,t — s) labels are omitted and J dfi[s, y \ t, x]stands 
for the standard' 16 ' path integral expression for the heat kernel. 

Having defined an Ito diffusion X(t) ,(5), (17) with the natural boundaries ,we are 
interested in the analogous (with respect to (35)) path measure Px 

P X [A] = j^dx J dnx[s, y\t,x]=J^ dfx x (36) 

The absolute continuity Px Pw implies the existence of the strictly positive Radon- 
Nikodym density, which we give in the Cameron-Martin-Girsanov form' 34,35 ! 



dpx 
djj, 



[s,y | t,x] = exp[f ^b(X(u))dX(u) - \ J* ±[b(X(u))] 2 du] (37) 



Notice that the standard normalisation appears ,if we set D — 1/2 which implies DA 
| A in the Fokker-Planck equation. 
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On account of our demand (9) and the Ito formula for Q(X(t),t) we have 



l - b{X{t))dX{t) = G(X(t),t) - Q(X(s), s) - f s du [d& + l -Vb]{X{u), u) 



(38) 



so that apparently 




[s, y\t,x] = exp[e(X(t),t) - Q{X{s), s)}exp [-^- f Q{X{u),u)du] 



J s 



(39) 



with = 2Dd t Q + DVb+ (l/2)6 2 introduced before in (10) , by means of the substitution 
of (4) in the Fokker-Planck equation. 

In case of natural boundaries at infinitythe connection with the Feynman-Kac for- 
mula (12) is obvious, and we have 



where the second integral refers to the path integration of the Radon-Nikodym density 
with respect to the conditional Wiener measure, see e.g. Refs.17,18,48-53,58. 

In the context of (40) and (12) we can safely assert that the pertinent processes 
(X(t) and W(t)) have coinciding sets of sample paths .The stochastic process "realizes" 
them merely (via sampling) with a probability distribution (frequency) different from this 
for the Wiener process W(t). 

The situation drastically changes ,if we wish to exploit the "likelihood ratio" for- 
mulas (37), (39) for diffusions confined between the unattainable (natural) boundaries, at 
least one of which is not at infinity. In view of the absolute continuity of Px with respect 
to Pi4/,we must be able to select a subset of Wiener paths which coincide with these ad- 
mitted by the process X(t), except on sets of measure zero (both with respect to Px and 



We face here a nontrivial problem of the existence of integral kernels for (Schrodinger) 
semigroup operators on a bounded or semi-bounded domain. The constrained version of 
the formulas (40), (12) should then integrate over a restricted set of Wiener paths:some 
of them must be totally excluded ,some must "avoid" certain areas (Wiener exclusion of 
Ref.42).This problem was solved in the context of the quantum statistical mechanics' 40-43,16 
where the Dirichlet boundary data for self- adjoint Hamiltonian are casually associated 
with the completely absorptive boundaries (cf.also Refs. 30-33, 36). The reason is clear ,if 
one lends priority to the Brownian motion proper, since there is no natural way for the 
standard Brownian motion (Wiener process) to prohibit it from reaching or passing any 
conceivable boundary, except for killing (stopping) the process, when it is going to happen. 

The most transparent way towards the localised Feynman-Kac representation of 
the Dirichlet (Schrodinger) semigroup on the (originally' 36 ' bounded) domain A is by 




(40) 



Pw). 



14 



introducing the first exit time T A for the Brownian path started inside A (a concrete 
sample path is labelled by u): 



T A (cu)=inf[t>0;X t (u;)eA} (41) 

Then ,the integral kernel of the (essentially self-adjoint on A, with the Dirichlet boundary 
data ) Hamiltonian H A = (—DA + (l/2_D)fi) A ,is to be given by the conditional expecta- 
tion: 

exp(-tH A )(s, y, t, x) = E y ^ s [exp[- J* Q(X u )du] ; t < T A | X t = X(t) = x] (42) 

which is an integration (40) restricted to these Brownian paths which while originating 
from y G A at time s,are conditioned to reach x G Aat time t, without crossing (but 
possibly touching) the boundary <9A of A. 

Another way to write down (42) is possible ,if we define a function a A (u) on the 
event set (event =sample path): 

= f 1, if X t (uj) G A for all * G [0, T] ^ 
1 0, otherwise 

Here a A is measurable with respect to the conditional Wiener measure, and then e.g. 
(35) can be replaced by the constrained path integral I 23 - 40 - 43 !; 

Pa[A] = j^dx J aA.{u))diJ, UJ [s,y \t,x] = J dfi A (44) 

where Ac A and oo is the sample path label (omitted in (35) to simplify notations). 

The analysis' 23 ' 40 ' 41 ' of special sets of Wiener measure zero is quite illuminating 
at this point. Namelythe integral (44) in addition to paths which are strictly interior 
to A, admits also paths which do touch the boundary <9A of A for at least one instant 
t G [0,T].Fortunately[ 23 ' 40 ' 41 l ,the P A measure of the set of such (unwanted) trajectories is 
equal zero. 

Let us now consider a diffusion X(t) , which lives in A and for which <9A is a natural 
boundary. Obviously no sample path of X(t) can reach (touch) d A in a finite time. By the 
absolute continuity Px <C Pw argument, we know that sets of Pw measure zero are the 
Px measure zero sets as well. Hence, (44) implies an apparent modification of (36): 

Px[A] = j^dx J a A (u)^^-[s,y \ t,x]diJi u [s,y,t,x] (45) 

applicable to the diffusion X(t) with the natural boundary dA. The Radon-Nikodym 
density is given by (39) and the path integral representation (42) is apparently valid for 
the involved (compare e.g. (4) again) Feynman-Kac kernel. 
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Except for the set of measure zero ,the process X(t) is characterised by the com- 
mon with the standard Brownian motion W(t) ensemble of sample trajectories,whose 
"realisations" by X(t) are appropriately (Cameron-Martin/Feynman-Kac) weighted. 

Although for each choice of the natural boundary <9A there is a unique diffusion 
, which respects it , we can devise a method of foliating the set of all considered diffusions 
into dynamically equivalent classes .We shall call diffusions dynamically equivalent, if they 
generate the same ,a priori given field of local accelerations bVb + DAb = Vf2 in their 
domain of definition.lt amounts to making a definite functional choice for V£l(x),x G 
i? x ,and then classifying all natural boundaries , which are consistent with this choice (let 
us emphasize that |6 2 + Wb = Q is to hold true modulo a constant). 

Following Ref.8 we can always consider a given ,unrestricted in R 1 Smoluchowski 
diffusion, as the reference one. Let Q(x) be its auxiliary potential, and Vf2 the induced 
field of local accelerations. What are the diffusions with natural boundaries ,which are 
dynamically equivalent to this reference one ? 

On purely technical grounds , the answer is simultaneously provided in the frame- 
work of Nelson's stochastic mechanics^ 4 ' 7 ' 8 ' 14 ' 18 ' 26-30 ' 48-53 ' 58 ] and of Zambrini's Euclidean 
quantum mechanics' 7,8 ' 50-53,58 ' 59 ' .The pertinent homogeneous diffusions belong to the over- 
lap of these two theoretical schemes and are uniquely specified by the nodal structure of 
stationary solutions of the Schrodinger-type equatiomwith D replacing H/2m in the orig- 
inal quantum evolution problem, and the Schrodinger potential being equal Q(x) modulo 
an additive (renormalisation) constant. The ground state process would correspond to the 
chosen Smoluchowski diffusion. 



Example: The notorious (albeit exceptional) harmonic attraction 
Let us consider the Sturm-Liouville problem on L 2 (R l ) 

2 2 

-DAip + = (46) 

The substitutions:^ 4 = uj 2 /AD 2 ,\ = e/uj,x = £/a give rise to the equivalent eigenvalue 
problem 

(~\\ + y)0 = ->4 (47) 

<Kt) = 1>(~) = I'M 
a 

with the well known solution (normalised relative to x) 

A n = n + - <-»• e n = (n + , n = 0, 1, 2, ... 
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M*) = UO = (^7^) 1/2e ^[- H n (0 (48) 

H = 1,H 1 = 2£,H 2 = 2(2£ 2 - 1) , H 3 = 4£(2£ 2 - 3) , ... 

Except for n = the solutions n (£) are not positive definite and change sign at nodes. We 
have 

n — , ipo(x) > , x G (— oo, +oo) 
n = 1 , ipi(x) > , x E (0, +oo) 
■01 (x) < , x G (-oo, 0) 
n = 2 , ^(x) > , x G (-oo, -l/v 7 ^) U (l/v 7 ^, +oo) 

ifa(x) <o,xe (-1/A +1/V2) 



n = 3 , 3 (x) > , x G (- y/3/2, 0) U (^3/2, oo) 

rl%{x) < , x G ( oo, -^3/2) U (0, ^3/2) 

and so on. It is convenient to continue further considerations with respect to the rescaled 
£ = ax variables, in view of the form — + y = H of the Hamiltonian predominantly 
used in the mathematical physics literature' 16 ' . To proceed in this notational convention 
it is enough to set x — * £ and D — * \ in the formulas (1)-(12) and thus utilize b = 
V0/0, n = \{b 2 + V6), Vfi = bVb + ±A6. 

Although in (1)-(12) we need G, 0* of the same sign, and p(x) to be strictly pos- 
itive, we can first make a formal identification 6 = 6* = (p n ,n = 0,1,2,... and notice 
that 

£ 2 1 

n = 0,6 = -£^fi = y- 2 
n = l,0i= £ — > s Zi = 

n = 2 , o 2 = — ^ t — > i 2 2 = 

2£ 2 -l s 2 2 

1 4£ £ 2 7 

n = 3 ' 63 = £ + 2e33^ 3 = y-2 

6bviously VfZ„ = £ for all n. Irrespective of the fact that each of b n , n > shows singular- 
ities,the auxiliary potentials are well defined for all x, and for different values of n they 
acquire an additive renormalisation — A n = — (n + |). 

The case of n = is a canonical' 16 ' 23 ! example of the Feynman-Kac integration, and 
the classic Mehler formula involves the Cameron-Martin-Girsanov density (39) as well. 
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Indeed^ 16 ' ,the integral kernel [exp(— Ht)](y, x) = k(y, 0, x, t) for H = — |A + (|x 2 — 
|) is known to be given by the formula: 

My, 0, t) = ^(1 - «-) - " 2 «*h^ " £^,1 W 

(e- H 'e)(i) = | 0, x, t)Q(y)dy 
where the integrability property 

/x — y 2 
k{y, 0, x, t)exp[ — - — }dy = 1 (50) 

is simply a statement (cf.(39)) pertaining to the transition density (4) of the homogeneous 
diffusion,which preserves the Gaussian distribution p(x) = (00*) (x) = -^exp(—£ 2 ). 

The case n=l automatically induces the (ergodic according to Ref.53) decomposition 
of the diffusion process into two independent, non-communicating components,each being 
confined between the natural boundaries (— oo, 0) and (0, oo), respectively. The pertinent 
processes have the same Feynman-Kac weight in the general expression (45) for their 
probability measures. Notice that we deal here with processes on the half-line, whose drift 
b± — | — £ has a singularity of the Bessel type,when the diffusion is to approach the 
point 0,see e.g. Ref.60 for a related discussion. It suggests that the construction of the 
probability measure on the half-line can be accomplished by directly starting from the 
Bessel process with natural boundaries at and oo. Namely ,the rescaled form of the 
backward Bessel generator 

L a = \^+ l -^V, (51) 

with x = v^C? a > 0, corresponds to the transition density of the diffusion with inaccessible 
boundaries 

pit; £ , = canst £ exp[-^M] J B (^) (52) 

where I a (ct) is the modified Bessel function. A particular choice of a = 1/2 i.e. I\/i{ot) = 
sinha provides us (in the notational convention (47)) with a conservative diffusion 
whose field of drifts is 6(£) = |.This diffusion can be directly compared (in the sense 

of Girsanov,see' 14 ' 18 ' 37 ' 38 ' 58 ' 60 ' 64 l to the unrestricted harmonic diffusion considered previ- 
ously:the drift transformation from | to | — £ induces a corresponding transformation 
of probability measures. Effectively it amounts to replacing in (45) the restricted (to the 
semi-axis) conditional Wiener measure by the conditional Bessel measure, which auto- 
matically respects the boundaries, and next evaluating the Radon-Nikodym derivative of 
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the harmonic measure with respect to the Bessel measure. The corresponding Cameron- 
Martin-Girsanov (likelihood ratio ) formula can be found in Ref.17. 

The decomposition into non-communicating diffusions with natural boundaries is 
characteristic of all n > solutions of (47). However all of them induce the same local field 
of accelerations Vfi(£) = £. 

Although the existence of the Feynman-Kac kernels (and thus of the transition den- 
sities and the diffusions themselves) is here granted ,it is generally not easy to give analytic 
expressions for them. However, in view of (44) ,the numerical simulation of each diffusion 
problem encountered before ,is definitely in the reach. Our discussion was basically one- 
dimensional and restricted to stationary cases, nevertheless extensions to time-dependent 
(non- stationary) processes and to higher dimensions (much of the outlined structure is 
preserved) are available. 
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